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Abstract 

Two "phase-shift equivalent" local NN potentials with different 
parametrizations, Reid93 and Nijmll, which were found to give nearly 
identical results for the triton by Friar et at, are shown to yield re- 
markably similar results for ^Li and ^^N in a (0 -I- 2)11^1 no-core space 
shell-model calculation. The results are compared with those for the 
widely used Hamada-Johnson hard-core and the original Reid soft- 
core potentials, which have larger deuteron D-state percentages. The 
strong correlation between the tensor strength and the nuclear bind- 
ing energy is confirmed. However, many nuclear-structure properties 
seem to be rather insensitive to the details of the NN potential and, 
therefore, cannot be used to test various NN potentials. 



1 Introduction 



Currently a large variety of nucleon-nucleon (NN) potentials exist. Despite 
the fact that all these potentials were designed to more or less reproduce 
the deuteron properties and the low-energy NN scattering data, they are 
quite different in details. In particular, the deuteron D-state percentage 
given by different potentials ranges from 4.4% for Bonn A |l| to 7.0% for 
Hamada- Johnson |p. Exact (or nearly exact) calculations for the A=3 and 
j4=4 systems |3|, |^ show that there is a strong correlation between the tensor 
strength of the one-boson-exchange NN potential and the binding energy. 
The NN potential which has a weaker tensor strength tends to yield a larger 
binding energy. 

In Ref.|^, however, it is found that when the NN potentials with different 
parametrizations are optimally fitted (with Xmin P^^ datum about 1) to the 
phase-shifts obtained in the "new and comprehensive (energy-dependent) 
partial-wave analysis", they result in nearly identical triton properties. It is 
pointed out that the tensor force is tightly constrained by the mixing param- 
eter ei, which can be determined very well in their multi-energy analysis. 
For example, three local NN potentials, Nijmll, Reid93 and AV18 with pa- 
rameters newly determined by their fit, have deuteron L'-state percentages 
of 5.64%, 5.70% and 5.65% respectively. 

It is the main purpose of this work to find out whether or not the above 
observation made in the study of the A=3 system also holds for larger nuclei. 
To this end, we will perform nuclear shell-model (SM) calculations for ^Li 
and ^^N using effective interactions derived from the Reid93 and Nijmll 
potentials. To set a scale for the comparison, the same calculations will be 
performed using the Hamada-Johnson hard-core potential (H-J) ||2[ and the 
original Reid-soft-core potential (RSC) both of which were widely used 
in nuclear-structure calculations in the past. We will adopt a no-core model 
space for which the effective interaction can reasonably be approximated 
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by the Brueckner G-matrix evaluated at appropriate starting energies. 
This avoids the difficulty encountered in doing the infinite summation over 
the core-polarization diagrams for the effective interaction in the presence 
of an inert core. Of course, for each nucleus, the effective interactions will 
be calculated in exactly the same manner so the difference reflected in the 
nuclear SM results is solely due to the difference in the NN potentials. 



2 Calculations 

The Brueckner G-matrix Q is calculated according to the following equation 
using the technique developed in Ref. |Q : 

G{lj) = Vi2 + Vi2 ^, . r^^l2, (1) 

LO - [hi + h2+ V12) 

where uj is the starting energy, V12 is the NN potential, h is the single- 
particle Hamiltonian with the one-body mean field V approximated by a 
shifted harmonic- oscillator (HO) potential u: 

/i = t + V~t-Fu-Vb = — + -mn'^r'^ - Vb- (2) 

2m 2 

The parameter Vq approximately represents the depth of the mean field V in 
which the two nucleons in the ladder-scattering processes are moving. This 
parameter can be absorbed into the starting energy by writing the G-matrix 
equation (||) as 

where uj' = lu + 2Vo and = 2^ + ^mQ'^r^. 

For the two-body G-matrix element {ab\G\cd) j^t , the shifted starting 
energy to' is taken as 

uj' = ea + eb + A, (4) 
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where e are the HO single-particle energies: 

2na + la + ^jhn. (5) 

Roughly speaking, the quantity A represents the interaction energy of the 
two nucleons bound by the HO potential. It is negative and depends on the 
two-nucleon state. In this work, A is regarded as an adjustable parameter 
whose value will be so chosen that we obtain a reasonable binding energy. 
It should be emphasized that A will only be adjusted when we go from one 
nucleus to another; its value is fixed when we consider the same nucleus but 
different NN potentials. We chose = 16 MeV, A = -38 MeV for A=6 
and nn = 14 MeV, A = -62 MeV for A=U. 

In the no-core model space, the starting-energy independent two-body 
effective interaction is simply the sum of ladder diagrams (G-matrix) plus 
the folded diagrams. When the starting energy for the G-matrix is prop- 
erly chosen, the folded diagrams have a small contribution to the effective 
interaction, which can, therefore, be well approximated by the G-matrix 
alone. 

Since we will be performing a (0 + 2)hi} calculation, the Pauli operator 
Q is defined as (A^i = 2ni -|- li and N2 = 2n2 + h are the principal quantum 
numbers for the two single-particle orbitals of the intermediate states in the 
ladder diagrams): 

For ^Li : Q = 1 for (iVi + N2) >4,Ni^ 0, and N2 ^ 0; 

= for (A^i + N2) < 4, A^i = 0, or N2 = 0. (6) 

For ^^N : Q = 1 for (A^i + N2) > 4, A^i > 1, and N2 > 1; 

= for (A^i + N2) < 4, iVi < 1, or N2 < 1. (7) 

Namely, in the ladder diagrams, the scattering into a two-particle state 
with ATi -I- iV2 < 4 is forbidden because these states will be included in the 
(0 -I- 2)?in calculation. In addition, for ^Li, the scattering into a two-particle 
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state with one particle in the Os orbit is also forbidden because the Os orbit 
is mostly occupied. For ""^^N, the scattering into a two-particle state with 
one particle in the Os or Op orbits is not allowed for the same reason. This is 
also shown in Fig.l. Notice that in Fig.l we have cut off the "wings" which 
extend out from the mostly occupied orbitals, at a large enough N, so as 
to obtain a converged result. It should be pointed out that our treatment 
of the Pauli operator Q is not exact; a more accurate definition is possible. 
Such a definition should be reasonable enough to allow us to do sensible 
calculations and meaningful comparisons. Here a more important issue is 
that the effective interactions for different potentials are calculated in exactly 
the same way. 

Once the G matrices are obtained, the Hamiltonian for the SM calcula- 
tion is written as: 

Hsu = (E - ^CM^ + E Gij + A (^HcM -l^^), (8) 

where Tcm is the center-of-mass (CM) kinetic energy and Hcm is the HO 
CM Hamiltonian. The parameter A is taken to be large enough (~ 10) to 
force the CM motion into its lowest configuration for all the low-lying states. 

3 Results and Discussions 

Our calculations are performed using the OXBASH SM code [^. The results 
for the ground-state energy, the low-lying positive-parity energy spectrum, 
and selected electromagnetic (EM) static and transitional properties in ^Li 
and ^^N are shown in Table I for the four NN potentials used in this work. 
The deuteron D-state percentage Pd for each potential is also listed. 

Before we proceed to discuss the results, it should be emphasized that 
our calculations involve some parameters, such as A and fiQ. The calcu- 
lated ground-state energy, which is quite sensitive to the choice of these 
parameters, should not be taken as an accurate result. Nevertheless, since 
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we perform exactly the same calculation for different potentials, the results 
are valid for the purpose of comparisons. 

In Ref.0], the Reid93 and Nijmll potentials have been shown to give 



almost the same results for the triton by solving the Faddeev equation |1C]. 
It is evident from Table I that the effective interactions calculated from these 
potentials also produce remarkably similar results for the more complex A=6 
and 14 systems. This is more obvious when we compare the results for these 
two potentials with those for the H-J and RSC potentials. Note that the 
similarity between the Reid93 and Nijmll results lies not only in the bulk 
properties, such as the binding energy, but also in more subtle observables, 
like excitation energies and EM transition rates. The latter reflects the 
equivalence in the wave functions for these two interactions. It is then safe 
to conclude that the Reid93 and Nijmll potentials are not only equivalent 
in describing the A=2 and 3 systems as shown in Ref. [^], they are also 
equivalent for the larger systems ^Li and ^^N studied in this work. 

For the A=3 and ^4=4 systems, it has been found that there is a strong 
correlation between the nuclear binding energy and the strength of the ten- 
sor component in the NN potential, characterized by the deuteron D-state 
percentage Pd 0, Apparently, the above correlation persists for the 
more complex nuclei ^Li and ^^N in a truncated-space SM calculation, as 
performed here. The calculated binding energies for the HJ and RSC po- 
tentials, which have larger deuteron Z)-state percentages, are much smaller 
than those for the Reid93 and Nijmll potentials. 

Note, however, that not all nuclear-structure properties are sensitive to 
the details of the NN potential. By examining Table I, one can see that 
the excitation energies for the low-lying states and some EM observables for 
different NN potentials are actually quite similar. This makes it difficult to 
test the NN potential in a nucleus by studying nuclear structure. 

On the other hand, besides the binding energy, there do exist a few 
observables which have a strong dependence on one or more particular com- 
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ponents of the NN potential. Two famous examples are the isospin splitting 
between the energies of the = 0^, T=0 and J'^ = 0^", T=l states in ^^O 
and the Gamow- Teller strength B{GT) between the ground state in ^^C and 
^^N. The former is sensitive to the tensor component in the NN potential 



[11 1 and the latter is sensitive to the spin-orbit and tensor components |T 
These have been studied in great detail in [13|. In Table I, one does see a 
sizable change in the B{GT) value from the ground state of ^^C, when going 
from one potential to the other. The B{M1) value is closely related to the 
-B(GT) value and also shows some change. However, both the -B(GT) and 
B{M1) values are not a monotonic function of Pd, suggesting that besides 
the tensor force strength, other differences between the four NN potentials 
also play a role. 

From Table I, we see that the calculated quadrupole moment for the 
ground state of ^Li also experiences quite significant changes when the NN 
potential is altered. In our truncated-space SM calculation with bare op- 
erators [i.e., e{p) = 1 and e(n) = 0], none of the NN potentials yields the 
experimental result. The smallness of the experimental quadrupole moment 
for ^Li and of the -B(GT) value from ^'^C seems to indicate that in a trun- 
cated SM calculation, one needs a weaker tensor force. This is consistent 



with the observations made in Ref.[13|, where some other quantities which 
are sensitive to the tensor force were considered. 



4 Conclusions 

In this work, we perform a no-core space (0 -|- 2)110. SM calculation for 
the ^=6 and 14 systems using effective interactions calculated from four 
different NN potentials. The Reid93 and Nijmll potentials, which were 
found to be equivalent in the description of the A=2 and ^4=3 systems, are 
shown to yield almost identical results, not only for the binding energy, but 
also for the energy spectrum and the EM static and transitional properties 
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for the larger nuclei considered in this work. It is verified that the NN 
potential with a weaker tensor component (i.e., smaller deuteron D-state 
percentage) tends to give a larger nuclear binding energy, consistent with 
the observations made in the study of the smaller ^=3 and 4 systems. Many 
nuclear-structure properties are found to be rather insensitive to the details 
of a NN potential and similar results are obtained for them, when different 
NN potentials are used. 

Among the various observables considered, there are a few which are 
sensitive to the various components of a NN potential, in particular, the 
tensor force. In order to reproduce the experimental results for the GT 
transition strength from the ground state of ^^C and the quadrupole moment 
of ^Li in a truncated space SM calculation, a weaker tensor force is needed. 
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Table I. The ground-state energies and excitation energies of low-lying states, in 
units of MeV, of ^Li and ^^N in a (0 -I- 2)hil calculation. The G-matrices are cal- 
culated in exactly the same way from different NN potentials including Hamada- 
Johnson (H-J), Reid-soft-core (RSC), new Reid-soft-core (Reid93), and new Ni- 
jmegen (Nijmll). For ^Li (^^N), a harmonic-oscillator basis with hfl = 16 MeV 
{hn = 14 MeV) is used. The GT strength S(GT) and Ml strength B{M1) (in 
fiff) is calculated for the transition 0^{1) 1^(0)- The magnetic dipole mo- 
ment fi (in hn) and electric quadrupole moment Q (in efm^) are calculated for 
the J^iT) = Ij^(O) state. Bare operators {gi{p) = 1, gi{n) = 0, gs{p) = 5.582, 
gs{n) = —3.826, e{p) = 1, e(n) = 0] are used in these calculations. In the Table, Pd 
is the deuteron Z)-state percentage. 
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Nucleus 


J'^{T) 


H-J 


RSC 


Rcid93 


Nijmll 


Expt. 




Pd 


7.00% 


6.47% 


5.70% 


5.64% 


N/A 




lf(0) 


-21.353 


-22.545 


-24.997 


-26.038 


-31.996 




3+(0) 


3.271 


2.946 


2.822 


2.840 


2.186 




0^(1) 


3.977 


3.691 


3.605 


3.566 


3.563 




2+(0) 


5.275 


5.038 


5.400 


5.448 


4.310 




2f(l) 


7.197 


6.816 


6.898 


6.906 


5.366 




1^(0) 


7.592 


7.617 


8.101 


8.205 


5.650 




2^(1) 


11.854 


11.530 


11.842 


11.957 






1^(1) 


12.840 


12.586 


13.065 


13.236 






0^(1) 


16.307 


16.345 


17.001 


17.204 






B(GT) 


5.327 


5.335 


5.387 


5.397 






B{U.l) 


15.21 


15.38 


15.53 


15.57 








0.840 


0.842 


0.848 


0.849 


0.822 




Q 


-0.265 


-0.208 


-0.168 


-0.163 


-0.11 




ii"(o) 


-56.306 


-75.944 


-87.562 


-87.087 


-104.659 




0^(1) 


2.176 


1.894 


1.702 


1.637 


2.313 




1 + {r\\ 


Z.59U 


1 Ai^ 


A /I 

Z.043 


AO 


A /I 

3.945 




2+(0) 


3.587 


3.374 


3.615 


3.594 


7.190 




2^(1) 


6.434 


6.660 


6.874 


6.801 


10.149 




3+(0) 


8.349 


9.227 


9.779 


9.757 


11.050 




1^(1) 


9.370 


9.492 


9.971 


9.971 


13.619 




B(GT) 


1.953 


1.243 


1.297 


1.305 







B{M1) 


5.380 


3.694 


3.823 


3.842 






M 


0.562 


0.495 


0.499 


0.499 






Q 


2.285 


2.155 


2.105 


2.109 
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Figure Caption 

Fig.l The Pauli operator Q for ^=6 and A=14 used in this work. 
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